Pion-Kaon correlations in central Au+Au collisions at y/s^N = 130 GeV 
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Pion-kaon correlation functions are constructed from central Au+Au STAR data taken at -^sjvjv = 
130 GeV by the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The results suggest 
that pions and kaons are not emitted at the same average space-time point. Space-momentum 
correlations, i.e. transverse flow, lead to a space-time emission asymmetry of pions and kaons that 
is consistent with the data. This result provides new independent evidence that the system created 
at RHIC undergoes a collective transverse expansion. 



Two-particle correlations for non-identical particles 
produced in heavy ion collisions are sensitive to differ- 
ences in the average emission time and position of the 
different particle species Such correlations in data 
taken at GANIL {^^^Xe+^8Ti at 45 MeV per nucleon) 
suggest delayed emission of deuterons with respect to 
protons 0. Correlation data from the SPS (Pb-Pb col- 
lisions at ^SNN — 17.3 GeV), and AGS (Au-Au colli- 
sions at y/SNN = 4.7 GeV) also suggest that the pion 



and proton average space-time emission points do not 
coincide 0, 0, 0; a partial explanation is that space- 
momentum correlations arise from the system's collec- 
tive expansion 0. For Au-|-Au coUisions at yj/vw = 130 
GeV, transverse mass spectra, elliptic flow, and deduced 
pion source radii suggest collective expansion in the 
transverse plane Such transverse flow may shift the 

average emission radii of different particle species by dif- 
ferent amounts. Also, different species may kinematically 
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decouple from the system at different times depending 
upon their interaction cross sections 0- In addition, the 
average emission time for a given species may be delayed 
significantly if produced dominantly through resonance 
decay. We construct pion-kaon correlation functions from 
Au+Au STAR data taken at ^/snn — 130 GeV and in- 
vestigate whether the pions and kaons arc emitted at the 
same average space-time position. 

Non-identical charged particles interact through 
Coulomb and strong interactions; for the pion-kaon case 
correlation effects are dominated by the Coulomb inter- 
action. To probe the tt-K separation, correlation func- 
tions C{k*) are constructed as the ratio of the k* dis- 
tribution constructed with particles from the same event 
(correlated distribution) divided by the k* distribution 
constructed with particles from different events (uncor- 
related distribution), k* is the magnitude of the three- 
momentum of either particle in the pair rest frame. 

For two particles initially moving towards each other 
the effects of the Coulomb and strong interactions are dif- 
ferent from those for two particles initially moving apart. 
The technique exploits this difference to study emission 
asymmetries ■ Pairs are divided into two groups, 

which represents either the case where the pions catch up 
with the kaons or the case where the pions move away 
from the kaons, depending upon the space-time separa- 
tion between pion and kaon emission points. Each sam- 
ple is used to construct two different correlation func- 
tions, C+(fc*) and C_(fc*), the sign index reflecting the 
sign of if-k*, with if the pair velocity and k* the pion 
momentum vector in the pair rest frame. If the average 
space-time emission points of pions and kaons coincide, 
both correlation functions are identical. If instead, pi- 
ons are emitted closer to the center of the source than 
kaons, pions with larger velocity will tend to catch up 
with kaons, and the Coulomb correlation strength will 
be enhanced compared to the case where pions are slower 
than kaons. Hence, the correlation function C+ will show 
a larger deviation from unity than C_ . Pairs can be sepa- 
rated according to the sign of k*-^^, ki^^^ and the k* 
projections onto three perpendicular axes in the longitu- 
dinally comoving system (LCMS) where the longitudinal 
component of the pair momentum vanishes j^]. The out 
axis parallels the pair velocity in the LCMS, the long 
axis is the beam axis and the side axis is perpendicular 
to the other two. r^^jj, r*^^^^, and rf^ng corre- 
sponding projections of the three- vector T**, the relative 
distance between the particle freeze-out points in the pair 
rest frame. Due to azimuthal symmetry and symmetry 
about mid-rapidity, (r^^J = {r^^J = 0. Thus C+/C_ 
defined with respect to the signs of fc*-^^ and k^^^^ must 
equal unity. If pions and kaons are not emitted at the 
same average radius in the transverse plane and/or at 
the same average time, C+/C_ defined with respect to 
the sign of k*^^ will deviate from unity, unless these two 



contributions cancel. Thus, one can probe the space-time 
separation between pion and kaon sources in the trans- 
verse plane. 

Charged particles are identified and tracked by the 
STAR Time Projection Chamber (TPC) H^- This anal- 
ysis selects the 12% most central collisions, i.e. the events 
with the largest multiplicity of particles. Selected parti- 
cles have pseudorapidity \ri\ < 0.5. The Au-|-Au colhsion 
point (primary vertex) is required to be within ±75 cm of 
the TPC mid-plane. The non-correlated pair background 
is constructed by mixing events whose primary vertices 
are also separated from each other by less than 10 cm. 

Pions and kaons are identified by measuring specific 
energy loss (dE/dx) in the TPC. When the momentum 
of pions and kaons exceeds 700 MeV/c, the dE/dx of 
both species becomes similar which compromises parti- 
cle identification. In addition, the pion and kaon samples 
are contaminated by electrons and positrons. The yield 
of each particle species in the momentum range where 
the energy losses coincide is interpolated {e~^ /e~' contam- 
ination) or extrapolated (kaon/pion separation) from the 
yields measured in the momentum range where there is 
good separation. In order to quantify the probability 
of correctly identifying a given species when the dE/dx 
bands overlap, four probabilities are calculated for each 
track: the chance that the particle is a 7r+ or tt", 
or , p or p, or or 6]. To be accepted as a 
pion or kaon the probability has to be > 60%. Tracks 
must point back to within 3 cm of the primary vertex; 
this removes a large number of secondary pions. Pions 
must have transverse momentum 80 MeV/c < pr < 250 
MeV/c and rapidity \y\ < 0.5, while kaons must have 400 
MeV/c <pT < 700 MeV/c, and \y\ < 0.5. 

Pion-kaon pair identification probability (product of 
both particle individual dE/dx probabilities) is required 
to be larger than 60%. Since the e~^-e~ pairs can distort 
TT~-K'^ and 'k~^-K~ correlation functions, the maximum 
probability allowed for a given pair to be e"'"-e~ is set 
at 1%, ensuring negligible contribution. Track pairs that 
share more than 10% of their TPC space points are dis- 
carded in order to avoid track merging errors. Two points 
are defined as shared if the probability of separating hits 
produced by them in the TPC is less than 99%. After 
selecting pion-kaon pairs, the correlation functions are 
constructed by taking the ratio of the k* distributions of 
pairs from the same event to the k* distributions of pairs 
from different events. 

Primary purity and momentum resolution effects are 
taken into account as described below. Primary purity 
is the percentage of primary pion-kaon pairs in all pion- 
kaon pairs satisfying all cuts. It is estimated to average 
77% for unlike sign pairs and 75% for like sign pairs. 
The lower limit for each is 54%. This number is the 
product of the probability of identifying both pions and 
kaons using the dE/dx information and the probability 
of excluding pions and kaons that do not originate from 
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points close to the collision vertex. Excluded pions in- 
clude decay products of strange hyperons and K^, and 
pions produced in the detector material. The fraction of 
secondary pions is estimated from the K^, A and pion 
yields in [ill . Detector simulations determine 

the relative reconstruction efficiency of pions from these 
different sources. Secondary kaons, being rare, are ne- 
glected. Assuming that the non-primary pion-kaon pairs 
are uncorrelated, the correlation function is corrected as 
Ctrue{k*) = {Cmeasured{k*)-l)/purity{k*) + l. The Sys- 
tematic error introduced by this correction is less than 
20%. 

The effect of momentum resolution depends upon the 
correlation function shape. Pion-kaon correlation func- 
tions are calculated from the pion and kaon momen- 
tum and space-time distributions, accountin g fo r both 
the Coulomb and strong interactions as in |l4|. The 
correlation function strength is calculated with the true 
momentum while the correlation function is binned as a 
function of k* smeared by momentum resolution. Mo- 
mentum resolution is estimated at the track level by de- 
tector simulations. The space-time distribution is cho- 
sen so that the main features of the measured correla- 
tion function are reproduced. The correction is obtained 
by comparing correlation functions calculated with and 
without momentum smearing. The correction enhances 
C{k*) by 20% (1%) for k* < 5 MeV/c {5 < k* < 10 
MeV/c), first and second bins in Figure 1, with a con- 
servative systematic error of ±100% on the correction of 
these two bins. 

The top panels of Figure 1 show the correlation func- 
tions for every combination of pion-kaon pairs. The 
agreement between unlike-sign {n^-K^ and n^-K^) and 
between like-sign {ir'^-K'^ and Tr~-K~) correlation func- 
tions is excellent. The middle and bottom panels show 
the ratios C+/C- for all pion-kaon pair combinations. 
C+ 1 C_ with respect to the sign of fc*-^^ and k*^^^ is unity 
within statistical errors in accordance with the require- 
ment that (r^rfj = (rf„„g) = 0. However, C+/C„ with 
respect to the sign of k*^,^^ is significantly larger than unity 
at small fc* when the interaction is attractive (t:~ -K^ 
and TT^-X") and significantly smaller than unity when 
the interaction is repulsive {j^^-K^ and Tr~-K~). These 
results indicate that pions and kaons are not emitted on 
average at the same radius and/or time. 

In order to understand the measured average space- 
time shift between pion and kaon sources, we compare the 
data with the RQMD (Relativistic Quantum Molecular 
Dynamic model and the Blast Wave Parametriza- 
tion (BWP) described in Ref. 0. BWP assumes that 
the system has undergone longitudinal and transverse ex- 
pansions, and provides the particle space-time and mo- 
mentum distributions at kinetic freeze-out. The param- 
eters, system outermost radius i? = 13 fm, freeze-out 
proper time r = 9 fm/c, emission duration At = fm/c, 
temperature T = 110 MeV, and transverse ffow rapid- 
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FIG. 1: Top panels: pion-kaon correlation functions C{k*), 
the average of C+(fc*) and C-{k*). Middle and bottom pan- 
els: ratio of the correlation functions C+ and C_ defined with 
the sign of the projections, fc*^^^ and ki^^g- Errors are 

statistical only. The horizontal axis of the ratios C+/C- for 
Kide {Kong) IS shifted by 1 MeV/c (2MeV/c) to separate the 
error bars. 



ity p{r) — 0.9{r/R) (with particle emission radius r) are 
consistent with fits to pion, kaon, proton and lambda 
transverse mass spectra and to pion source radii Q . The 
hadronic cascade model, RQMD, also generates trans- 
verse flow through rescattering of hadrons Indeed, 
turning off hadronic rescattering within this model shuts 
off transverse flow 16] . In addition, RQMD includes 
contributions from resonance decay, such as lu, rj and 4>, 
which shift pion and kaon emission times. 

Figure 2 shows correlation functions C{k*) and ratios 
C+/C- measured, from BWP, and from RQMD with 
and without hadronic rescattering. The calculated cor- 
relation functions use model space-time and momentum 
distributions as described in 0|, with pion and kaon 
kinematic cuts chosen to match the data. The correla- 
tion function are calculated for like-sign and unlike-sign 
pairs The small wiggles in the calculated C+/C_ ratios 
for k* < 20 MeV/c arrise from statistical uncertainties. 
RQMD and BWP are in qualitative agreement with the 
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FIG. 2: Comparison of the correlation functions between 
data and model. Upper panel, C{k*) correlation function. 
Lower panel, ratio C+/C- with respect to the sign of k^ut- 



measured correlation functions. Turning off rescatter- 
ing in RQMD leads to a strong correlation, which im- 
plies that the pion and kaon sources are too small. On 
the other hand, RQMD reproduces qualitatively the ratio 

The effect of source size and source shift is disentan- 
gled by simultaneously fitting the correlation functions 
C+ and C-. In order to insure that the detector ac- 
ceptance is matched, the particle momenta are taken 
from experimental pion-kaon pairs constructed by mixing 
events that pass all the cuts. The particle positions are 
set such that the distribution of the relative space-time 
separation between pions and kaons in the pair rest frame 
is a three dimensional Gaussian. The free parameters 
are the Gaussian mean, (Ar*^^) = (r*„t(7r) - r^^tiK)) 
((Ar*-^J = {^r*i^ng) = 0) and the Gaussian width, 
a — <Tr' — CTr' = cr^' . Both fit parameters from all 

oxLt side long 

four correlation functions are in agreement within statis- 
tical errors; combined they give a — 12.5±0.4t|2 fmand 
(Ar*„4) = -5.6±0.6tli with a x V dof = 134.5/110. 
Systematic errors are estimated from the discrepancy be- 
tween the four correlation functions, the dependence on 
the input momentum distribution, the uncertainties on 
primary purity and the fit range dependence. This -5.6 
fm in the pair rest frame becomes in the lab frame -3.9 
fm (5.4 fm/c) if emission difference is space (time) only. 

The parameters a and (Ar*^^) may be extracted 
directly from BWP or RQMD by constructing the 

r* = ^Jir*out)^ + iK^dey + i^tongy ^ud r*„, distribu- 
tions. However, neither RQMD nor BWP T** distribu- 



Data 12.5±0.4j:^'2 ~5.6±0.6t\l 134.5/110 
RQMD 11.8 ±0.4 -8.0 ±0.6 205/54 

RQMD 5.8 ±0.1 -2.0 ±0.3 

no rescattering 

BWP 9.9 ±0.1 -6.9 ±0.3 1020/118 



940/54 



TABLE I: Fit results using a three dimensional Gaussian dis- 
tribution in the pair rest frame. For the data, the first error is 
statistical and the second systematic. The errors on the model 
calculations are calculated by rescaling the x^ distribution by 
the minimum value of x^/dof. 



tion is close to a three dimensional Gaussian. Thus, to 
compare models and data fairly, the correlation functions 
calculated from RQMD and BWP are fitted in exactly 
the same way as the data. The extracted fit parame- 
ters are compared to the data in Table 1. The large x^/ 
dof values arise because the tails of the 'r'* distributions 
of RQMD and BWP are not well-described by a three 
dimensional Gaussian in the pair rest frame. The data 
appear to be insensitive to these tails due to larger sta- 
tistical errors. 

Consider BWP. At an emission point, the fluid velocity 
(increases with radius) and the thermal velocity (com- 
mon freeze-out temperature T for all species in fluid rest 
frame) combine to give the observed particle velocity \^ . 
If the source does not expand, the relative emission prob- 
ability for given will track the fireball spatial density. 
If the source expands but T = 0, particles with 1? will 
come from the single point where the fluid moves with l7. 
At T 7^ and for constant density and unlimited fireball 
size, the spread of thermal velocity smears this emission 
point to a nearly spherical volume whose size increases in- 
versely with particle mass. This volume must be folded 
with a realistic fireball spatial density distribution, re- 
moving contributions from large radial distances. Thus, 
effective centers of emission regions are shifted towards 
smaller radii. For our rrit/T [rrit = transverse mass, and 
TOt cx m at given V), the relative shift of pions and kaons 
is small 2] but measurable. There is also an emission 
time separation: BWP has kinetic freeze-out at fixed lon- 
gitudinal proper time r = ^/t"^ — z^, so the larger size of 
effective pion source yields emission at later laboratory 
times t. Thus pions are on average emitted closer to 
source center and later in time than kaons. 

In the RQMD model, pion and kaon sources are 
also spatially shifted when transverse flow builds up by 
hadronic rescattering. Even when rescattering is turned 
off, resonance decays delay the pion average emission 
time and increase the apparent size of the source. 

Our results show that pions and kaons are not emit- 
ted at the same average space-time position for Au+Au 
collisions at ^/smn = 130 GeV. The data are consistent 
with BWP and RQMD, i.e. with a system whose domi- 
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nant feature is a transverse collective expansion. These 
results significantly challenge models that attempt to ex- 
plain pion, kaon and proton spectra by purely initial state 
effects 0, 01 . Such an analysis may also be used to 
probe at what transverse momentum soft processes (ex- 
panding system) give way to hard processes since the 
space-time emission pattern will substantially change at 
that momentum. 
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